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Abstract 

We present a new perturbative approach to QCD based on the use 
of quark composites as fundamental variables. The composites with 
the quantum numbers of the nucleons are assumed as new integration 
variables in the Berezin integral which defines the partition function, 
while the composites with the quantum numbers of the chiral mesons 
are replaced by auxiliary bosonic fields. The action is modified by the 
addition of irrelevant operators which provide the kinetic terms for 
these composites, and the quark action is treated as a perturbation. 
The resulting expansion has the quark confinement built in. As a first 
application we investigate the pion-nucleon interaction. 
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1 Introduction 



We can divide the calculations in QCD in two classes: Those which aim at 
clarifying the confining mechanism, and those devoted to the description of 
hadronic physics. One of the present the difficulties in this field stems from 
the fact that the former calculations are done in a framework where it would 
be very awkward, if not impossible, to perform the latters. More generally, 
the actual understanding of the theory requires a nonperturbative approach 
for low energy and perturbative methods for high energy, which until now it 
has been impossible to unify in a unique scheme. 

The strategy we adopt to overcome the impasse is to use quark composites 
with hadronic quantum numbers as fundamental variables An earlier 
attempt in this spirit can be found in but it is restricted to the strong 
coupling. We avoid such limitation by means of new technical tools. 

There is a fundamental difference between trilinear and bilinear compos- 
ites. Surprisingly enough some trilinear composites, in particular the nucleon 
fields, can easily be assumed as integration variables in the Berezin integral 
which defines the partition function, since they satisfy the same integration 
rule as the quark fields. As a consequence the free action of these composites 
is the Dirac action ||. 

Also the mesonic composites can be introduced as integration variables, 
but the resulting integral does not reduce in general to a Berezin or to an 
ordinary integral, and it is very complicated. In particular the propagator is 
not the inverse of the wave operator. It has then required some effort to find 
out an operator containing the free action of the chiral composites M, but 
once it has been constructed, we can circumvent the difficulty of the integral 
over the chiral composites by replacing them by auxiliary fields by means of 
the Stratonovich-Hubbard representation ||. 

Since the free composite actions are irrelevant operators, we can freely add 
them to the standard action, and derive a perturbative expansion by treating 
the quark action as a perturbation. The expansion parameters are the inverse 
of two dimensionless constants entering the definition of the nucleon and 
chiral composites, and of the number of components of the up and down 
quarks. This number, called Q, is equal to 24. 

Our approach has two quite desirable features. Firstly it is compatible 
with a perturbative as well as nonperturbative regime of the gluons wr to 
the gauge coupling constant. For this reason we adopt a regularization on 
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a Euclidean lattice, which is the only one suitable for both cases. This 
choice seems also natural dealing with composites. The other remarkable 
aspect is that the quark confinement is built in. This is due to the fact that 
our perturbative expansion is of weak coupling for the hadrons, but, as a 
consequence of the spontaneous breaking of the chiral symmetry, of strong 
coupling for the quarks. To any finite order, the quarks can then move 
only by a finite number of lattice spacings, and therefore they can never be 
produced in the continuum limit, at variance with the standard perturbation 
theory, where the the quarks are propagating particles at the perturbative 
level, and the confinement is essentially nonperturbative. Let us emphasize 
that the hopping character of the quark propagator should not generate any 
confusion with a strong coupling expansion in terms of the gauge coupling: 
As stated above, the gauge coupling constant can have any value required by 
the dynamics, while mesons and nucleons do move in the true continuum. 

One can wonder whether the Wilson term to avoid spurious quark states 
is at all necessary in our approach: Since the quarks have no poles whatsoever 
to finite order, why should we worry about the spurious ones? In any case we 
can safely assume the Wilson term for the quarks of subleading order in our 
expansion parameters, so that we can ignore it in a first order calculation. 

In the previous works the nucleons and the chiral mesons were discussed 
separately. Here we unify their treatment including the em interactions, so 
as to be able to evaluate strong corrections to em processes. We derive a 
perturbative expansion in terms of the nucleon composites and the auxiliary 
fields which replace the chiral composites. The numerical coefficients of this 
series are given in terms of integrals over the gluon and quark fields. As 
an application we investigate the nucleon-pion interaction, getting to lowest 
order a contribution to be the one-pion exchange of the Yukawa model. 

The paper is organized as follows. In Section 2 we define the modified 
partition function we will use. In Section 3 we define the hadron composites 
and their actions. In Sections 4 and 5 we report for the convenience of the 
reader the results, obtained in the previous works, concerning the introduc- 
tion the nucleon fields as integration variables, and of the auxiliary fields for 
the chiral composites. In Section 6 we derive our perturbative expansion. 
In Section 7 we show that the quarks do not propagate at the perturbative 
level. In Section 8 we evaluate the nucleon-pion interaction and in Section 9 
we conclude our paper. 
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2 The modified partition function 



We assume the modified partition function 

Z = J [dU] J [dA] exp[-S YM - S M ] J [dXdX] exp[-S N - S c - S Q ], (2.1) 

where Sym and Sm are the Yang-Mills and Maxwell actions, Sq is the action 
of the quark fields and Sn and Sc are irrelevant operators which provide 
the kinetic terms for the quark composites with the quantum numbers of 
the nucleons and the chiral mesons. A is the quark field, the gluon field is 
associated to the link variables and A^ is the photon field. The reason 
why we have introduced the photon field explicitely, is that it is associated 
to different link variables in the interactions with the quarks and the com- 
posites. Finally differentials in square brackets are understood, as usual, as 
the product of differentials over the sites and the intrinsic indices, and [dU] 
is the Haar measure. 

All the elementary fields live in an anysotropic euclidean lattice of spacing 
a in the spatial directions and a t in the time direction, whose sites are identi- 
fied by fourvectors x of spatial components x^ = 0, ...N and time component 
t = 0,...2N t , and satisfy periodic boundary conditions, with the exception 
of the quark fields which are antiperiodic in time: X(x) = X(x + Nek) = 
—X(x + 2N t et), ek and et being the unit vectors in the k and time directions. 

The simplest way of unifying the treatment of nucleons and mesons is to 
put them at different, complementary sites of the latttice. A possible op- 
tion, which respects the hypercubic symmetry, is to put them on the sites 
of two dual sublattices, defining the parallel transport along the major di- 
agonals [0]. Here we make a different choice, possibly more suitable for the 
hadron termodynamics. We restrict the nucleons/mesons to the odd/even 
sites in one direction ( which is natural to assume as the time direction), 
while letting them to occupy the whole remaining three dimensional volume. 
The hypercubic symmetry is lost, but it remains valid at least for the free 
composites, if we assume the temporal spacing half the spatial one. We will 
keep however the spacings unrelated having the hadron thermodynamics in 
mind. So the nucleon and chiral composites are defined respectively at the 
sites n = (x fc , 2t + 1) and c = (x^, 2t),t — 0, ...N t , which identify the nucleon 
and chiral lattices. Accordingly we introduce specific notations for the sums 
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over the whole lattice or its sublattices 
(f,g) = a 3 a t Y<f( x )g( x )i 

X 

(f,g) n = a 3 2a t J2f(n)g(n), (f,g) c = a 3 2a t J2f(c)g(c). (2.2) 

n c 

The quark fields \® a with color, isospin and Dirac indices a, r and a 
respectively, are related to the up and down quarks according to 

AL=<, A^ = C (2.3) 

Their action is 

S Q = (\,Q\), (2.4) 

where 

Qx, y = M Q 5 xy - — ^2(r q - -ik)U k (x)v k (x)5 x+ek) y 
Aa k 

-7^-^2( r q-let)U e t(x)v et (x)5 x+eetj y, e = ±l, (2.5) 



with 

M Q = m q + r q (- + — . (2.6) 
We have adopted the standard conventions 



k G {-3, ...,3}, e_ fe = -e fc , e_ t = -e t , 
7-fc = -7fc> 7-t = -7t> 

[/_*(*) = C/+(x - e k ), U- t (x) = U t + (x - et), (2.7) 

and the corresponding ones for the link variables associated to the photon 
field 

v kr (x) = exp[iaq T A k (x)}, v tT {x) = exp[ia t q T A t (x)]. (2.8) 
q T is the quark electric charge 

2 1 

qi = -e, q2 = — 776, e = electron charge. (2.9) 
o o 

The spurious quark states are prevented by the Wilson term with Wilson 
parameter r q . 
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Since we will assume the nucleoli and chiral composites as fundamental 
variables, we must separate the quark action into the terms which live only 
in the nucleon and chiral lattices and those which couple them 



Sq — SQ nc + SQ cn + Sq u + Sq c (2.10) 



where 



Sqtic — 


a 3 a t ^2 \(n)Q n , c \(c 

nc 


Sqcu = 


a 3 a t ^2\{c)Q c , n \(n 

cn 


$Qn = 


-(A, Q\) n , 


Sq c = 


^(A,QA) C . 



(2.11) 
We then have 

Z = J [dA] exp[-S M ] J [c/Ac/A] n exp(-S' A r - S Qn ) J [dU] exp[-S Y M] 

J [d\d\] c exp [-S c - S Qc - S Qnc - S Qcn ] , (2.12) 

where 

[dXdX] = [dXdX] n [dXdX] c , (2.13) 

The formalism for a nonlinear change of variables in Berezin integrals has 
been developed in refs. PJ, 0. In the following section we will report what 
is necessary to arrive at the formulation of our perturbative expansion. 



3 The hadron composites and their actions 

The nucleon composites are 

2 1/2 

Ipra = a e a 1 a 2 a a S TT2 € TlT3 ('y5-)' fl ) aai {Cjn)a2a 3 K\ai K%a 2 K$a 3 - (3-1) 

In the above equation and in the sequel the summation over repeated 
indices is understood, C is the charge conjugation matrix, and ipi ai ip2a are 
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the proton, neutron fields. It is easy to check that with the above definition 
they transform like the quarks under isospin, chiral and 0(4) transformations. 

We should explain why we have included in the definition the cubic power 
of the lattice spacing a and the parameter k]y. The cubic power of a param- 
eter with the dimension of a length, say / 3 , is necessary to give the nucleon 
field the canonical dimension of a fermion field. At the same time a power 
of the lattice spacing at least cubic is necessary to make the kinetic term 
(with the Dirac action), irrelevant. We have written for later convenience I 3 
in the form k]( 2 a 3 , where fcjv is dimensionless and it must not diverge in the 
continuum limit. 

The nucleon free action is the Dirac action 

S N =$,Nj>) n , (3.2) 

whose wave operator N is given by 

/3 1 \1 1 

m N + T N y- + —J 5 nun2 - — J2( r N ~ 7fe)^ / fc(^l)^n 1 +e fc ,n 2 



N 

1 * 111, 112 



Y,( r * - 7 e t)V r rt(ni)£ ni+ 2e etl „ a . ( 3 - 3 ) 

In the above equation is the Wilson parameter 

< r N < 1, (3.4) 

Ttijsj is the mass of the nucleon and V M is the link variable associated to the 
e.m. field acting on the nucleons 

Vfc T (n) = exp(ie T aA k (n)), 

V tT {n) = exp(ie T a t A t (n)) exp(ie T a t A t (n + e t )), e x = e, e 2 = 0, (3.5) 

with the standard conventions. Needless to say, the coupling with the em 
field is necessary to preserve the abelian gauge invariance of the QCD action. 
The chiral composites are the pions and the sigma 

7r± = i k w a 2 \ r y 5 —z=(r 1 ± ir 2 )\, tt q = i k n a 2 A7 5 r 3 A, a = k n a 2 \\, (3.6) 
V2 

where 75 is assumed hermitean, the r^'s are the Pauli matrices and the factor 
a 2 k n has been introduced with the same criterium used for the nucleons, so 
that also k n must not diverge in the continuum limit. 
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The chiral transformations over the quarks 



i 



5\ = -7 5 r-5A (3.7) 
6X = ~A 75 f-a (3.8) 

induce O (^-transformations over the mesons 

8a = a ■ 7T (3.9) 
Sjf = -aa. (3.10) 

For any real vectors tt,x, we adopt the convenction 

7T • X = K+X- + 7T-X+ + TTOXO- (3.11) 

Since for massless quarks the QCD action is chirally invariant, the ac- 
tion of the chiral mesons must be, apart from a linear breaking term, 0(4) 
invariant. It must then have the form 



S, 



c 



i(vr, Ctx) c + -(a, Ca) c - — (Vttm, a) 

2 2 n z v / < 



(3.12) 



The factor yVt has been introduced for later convenience. The heuristic 
considerations which led to the choice of the wave operator C are based on 
experience with simple, solvable models, and can be found in [|]] 

C ( w ) = a ' 2 2 -nf 2 - ( 3 - 13 ) 
a z L> z — ft 

C is given in terms of the covariant laplacian T> 2 whose action on a function 
g is 

(V 2 g) (c) = ^ Y.\Mc)g{c + e k ) - g{c)} + £K(c)#( C + 2e ei ) - 0(c)] 

(3.14) 

w being the link variable associated to the em field acting on the charged 
pions 

w k (c)Tc± = exp(TieaA k (c))n ± , 

w t (c)-K± = exp(Tiea t A t (c)) exp(Tiea t A t (c + e t j)ir ± . (3.15) 
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Notice that the 7r± has negative/positive charge. A posteriori the choice of 
the wave operator can be justified by observing that in the Stratonovitch- 
Hubbard representation the kinetic term of the auxiliary fields is related to 
C^ 1 , and therefore the present form generates in the simplest way a Klein- 
Gordon operator. Let us stress that if the 7r' and the a were ordinary bosons, 
rather than even Grassmann variables, with the above action they would be 
static. This is not a paradox, because the propagator of even Grassmann 
fields is not the inverse of the wave operator. 

The irrelevance of the nucleon action is ensured by the insertion of the 
appropriate powers of the lattice spacing in the definition of the nucleon com- 
posites. For the chiral action instead the request of irrelevance constrains also 
the parameter p which cannot vanish in the continuum limit. The depen- 
dence on a assumed below for the breaking parameter to, ensures also the 
irrelevance of the chiral symmetry breaking term. 



4 The nucleon composites as integration vari- 
ables 

In this Section we report the results we will need on the change of variables 
for trilinear composites. 

There are 8 nucleon field components at any site n. Let us denote by 
^(n) their product in the order 

= ipu(n)ip 12 (n)...ip23(n)^24(n). (4.1) 

These fields obey the same multiplication rules of the quark fields. We can 
then define for them an integral 

/ Y[dip T ,a{n)^(n) = 1, all the other integrals vanishing, (4.2) 

which obeys the same rule as the Berezin integral for the quarks 

/ { [ d\ c ! ra (n)A(n) = 1, all the other integrals vanishing. (4.3) 

a,r,a 

In the latter equation A(n) is the product of all the quark components at the 
site n in the order 

A(n) = P(A 11 (n))...P(A 14 (n))P(A 21 (n))...P(A 24 (n)), (4.4) 
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with 

P(X Ta (n)) = \ l Ta {n)\l a {n)\Un). (4.5) 

It is then easy to check that for an arbitrary function which depends on the 
quark fields only through the ip, we have, with the appropriate ordering for 
the differentials 

/ II d\%{n)g{^{\{n))) = J f l[d^ a (n)g^(n)), (4.6) 

a,r,a T,a 

provided J, defined by 

*(n) = JA(n), (4.7) 

is different from zero. This is in fact true for the nucleons and, for the color 
group SU (3) in 4 dimensions 

J = k% a 24 2 22 • 3 3 • 5. (4.8) 

This result justifies the choice of the Dirac action as the free action for 
the nucleons. Let us in fact define the free nucleon partition function 

Z N = J[dXX] n exp[-S N ]. (4.9) 

Notice that the integral is over the quark fields in the nucleon lattice only. 
Let us now consider the nucleon-nucleon correlation functions 

< i> T MMy) >= it I [d\\}n^ Ta (x)^(y) exp[-SW]. (4.10) 

Zj N J 

Since the integrand is a function of the quark fields only through the ip, we 
can assume them as integration variables getting 

< ^ Ta (x)^{y) >= —J^y^^aiy - x), (4.11) 

namely the free correlation function of two Dirac particles. It should be noted 
that this result does depend neither on the value of the jacobian, provided 
it is different from zero, nor on the value oi k^. Inverse powers of this latter 
constant, however, will appear in the perturbative expansion of QCD, due to 
the presence of Sq in the total partition function. 
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It is remarkable, and it is at the basis of our perturbation theory, that 
more general formulae hold, which involve the quark as well as the nucleon 
fields. In fact the following substitution rule holds in Berezin integrals 



VT»-V>7 ro (n) ^h(n)-ipi m ,{n), (4.12) 

where the symbol ~ does not refer to any approximation, but it means that 
the equality holds only under the Berezin integral, and the /' are numerical 
coefficients called transformation functions. We have used the shorthand 
notation hi = (a i; Tj, clj), hi = (ai,Ti,azi), I = (r,a),I = (r, a). Restricting 
ourselves to integrals trilinear in the quark and antiquark fields 

II d \r,a {n) d\ a ^ a {n)g(ip (n) , ip (n) )\ i {n)\ 2 {n)\ 3 (n) \ hl (n) X h% (n) \ hs (n) 
flh 3 h2hJihih2h 3 J 2 / nd$ Tta (n)d^ T)a (n)g($,il>)^n) ipi(n). (4.13) 



The transformation functions in such a case are 

(4.14) 

where 



1 



3 7.-1/2 



h n~ 6 h 

9o 



U2tt:i 



0103 • (4.15) 

The functions h are totally symmetric wr to the exchange of any pair of num- 
bered indices. The appearance of /%^ 2 in the rhs shows how an expansion 
in this parameter is generated. 

It should now be clear in which sense we can talk of a change of variables. 
Even though the quarks cannot be expressed in terms of the composites, we 
only need to invert monomials in the quark and antiquark fields of order 
3-integer, and this can be done according to Eq.( |4.12| ). 
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5 The auxiliary fields for the chiral compos- 
ites 



In this Section we review the properties of the chiral action necessary for the 
further developments. To this end we define the chiral partition function Z c 



Z c = J [dA] exp[-S M ] J [d\d\) c exp[-S c } 



(5.1) 



where there appear the chiral composites and the Maxwell field only. Also the 
chiral composites can be assumed as integration variables, but the resulting 
integral is impracticable. Since Sc is quadratic and the wave operator C 
is negative definite, the difficulty can be overcome by introducing auxiliary 
fields by means of the Stratonovich-Hubbard transformation 



exp[-Sy 



(det C(A)) 
1 



-i 



dx 



cxp 
exp 



(£(a 4 C)~ 1 x) c +(0, (a 4 C)"V) 



a 2 a 2 v /( 



(5.2) 



In the above equation we have ignored, as we will do in the sequel, field 
independent factors. Now the chiral partition function 



Z c = J[dA]exp[-S M -Trln^O)- 1 ^))] 
dx 



I 



dd> 



exp <; -p 



(£ (a 4 C)- 1 x) c + (0, (a 4 CTV) ( 



[d\d\] c exp 



-y y 

—P 2 {X, tt)c + — (pV + v^m, o) 



(5.3) 



describes bosonic fields interacting with the quarks. To get the effective 
action of the chiral mesons we integrate out the quark fields with the result 



Z c = J[dA]exp[-S M - TrMCiOr'CiA))] J 



dx 



exp[-S x ], 
(5.4) 
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where 



-| J> {a 2 *; 2 [(v^m + p 2 0(c)) 2 + p 4 x{cf] } • (5-5) 



Since Q is a rather large number we can apply the saddle-point method and 
evaluate the partition function as a series in inverse powers of this parameter. 
The minimum of S x is achieved for x — and 



= 



y/a2a t p 



where 



£ = —y/a2a t m. 
2p 



The quadratic part of S x is 



where 



0-0 



is the fluctuation of the field 6 and 



m = 



m , 



a 2 6 ' 

2p 2 v^+F 



with 

Therefore, if we assume 



m 



-p^ t aml + ° {1/VUh 



(5.6) 
(5.7) 

(5.8) 
(5.9) 



(5.10) 
(5.11) 

(5.12) 
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in the continuum limit 



m_ 



m. 



2p 2 



(5.13) 



and the propagator of the pion field to leading order in 1/yVL turns out to 
be the canonical one 



< 7r/ l (ci)7r fc (c 2 ) >= 



1 



1 



a 3 2a t \~V 2 



mi 



(5.14) 



Cl,C 2 



while the a is unphysical because its mass is divergent. Notice that, in 
analogy to the case of the nucleons, the constant k n is uninfluent at this 
stage, but its inverse powers will appear in our perturbative expansion of 
QCD. 

The nucleon action Sn gives an exactly free propagator. Instead, the 
chiral action contains a residual interaction which can be obtained by 
expanding the In in inverse powers of Q. Here we report the first terms, 
given in [Q (a factor — is missing in front of and appearing in 
Eq.(63) of this paper) 



'2at 
a 



H? 2 + 3(x) 2 ) c 



+ 



4Q \b 



p\ 4 2a t r 



\ d\ + (X 2 , (X 2 - 6 9 2 )) c ] + 0(fi- 3 / 2 ).(5.15) 



6 The perturbative expansion 

In this section we derive a perturbative expansion for the partition function 
of QCD in terms of the nucleon composites and the auxiliary fields. We 
start by the Stratonovich-Hubbard transformation to get rid of the chiral 
composites 



exp[-Sc - S Qc ] 



(detC(A)r 1 
A exp 



dx 



^7T 



exp[-S A 



-(\,(D + Q H )\) C 
a 



(6.1) 
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where 



A 
D 



(ndetD,)- 1 



ak n [, 



p'(p + VQm + ip 2 7 5 r ■ x 



(Q 



7 E( r - 7fc)^fc(ci)wfc(ci)5 Cl+efeiC2 . 
4 fc 



(6.2) 



We have absorbed the term Mq in the parameter m which now reads 

l r, 1 



m = 



-Mr 



p V 2a < "" _7r '-Vl!/,- 

The partition function is so transformed into 

[dA] exp[-S M - Tr \n(C-\0)C(A))} / [dAdA] n exp[-S 



(6.3) 



'AT 



dx 



de 



Zexp[— S 



(6.4) 



where 



1 



[dU] exp[-S Y M] exp[-S Qn ] A / [d\d\] 



exp 



-(A, (D + Qh)X) c - S Qnc - 

(X 



(6.5) 



To perform the integral over the quark fields we expand X with respect to 
Sq. We note that only terms with an equal number of factors SQ nc and Sqcu 
can contribute. Moreover, after the integration in the chiral lattice there 
remain in I only polynomials of the quark fields in the nucleon lattice. They 
will contribute, according to Eq. (|4.12|) , only if they are of order 3-integer 
both in A and A. Therefore 

r 00 1 r - 

1 ~ J[dU]exp[-S YM }Y: E (- 1 ) S 7^y(^) SA J /[^] C 



{SQnc) r (S Qcn y exp 



i=o r+s=3i 

-(\,{D + Q H )X) ( 
a 



E E ^s. (6.6) 

i=o r+s=3i 
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We remaind the reader that the symbol ~ does not refer to any approxima- 
tion, but it means that the equality holds only under the Berezin integral. 
The functions I r>s can be evaluated by expanding the exponential wr to Qh 



t=0 



Since only even powers of Qh can contribute 



(6.7) 



J, 



v,/,u 



[dU}ex V [-S YM }(-iy 



(SQncY(SQ c 



(r!) 2 s!(2t)! 

1 _ \ 2t 
-(X,Q H X) C ) exp 
a / 



(s Qn y a / [d\d\] 



-(A, DX) C 



(6.8) 



Finally assuming the nucleon composites as integration variables we get 
the desired perturbative expansion of the partition function of QCD in terms 
of the nucleon and chiral fields 



Z = j '[dA] exp[-S M -Tr In (C^^C '(A))] j '{d^dvfj}, 



dx 
s/2n 



'2tt 



l r , s exp[-S N - S x ]. 



(6.9) 



c i=o r+s=3i 



The integrals over the quark and gluon fields are relegated in the functions 
X r s . Note that we do not need to treat the gluon field perturbatively: If and 
where this can possibly be done remains here an open question. 

At this point the role of fc^ 1 , k~ l as expansion parameters should be clear. 
Since the integral over the quarks in the chiral lattice is proportional to D^ 1 , 
namely to fc~^ 2 * +r - ) , and because of the dependence of the transformation 
functions (gig ) on k N , X r ^ t ~ k'/^'^k-^^. 

We conclude this Section by the evaluation of the correlation functions. 
To this aim we observe that the Stratonovich-Hubbard transformation can 
be written in the form 



7T hl {C l )...'K hs {C s 

exp[— S x }(a2atp 2 



exp[-Sc 




S Qc ] = (detC(A)y 1 J 



dx 
s/2tt 



d 



dXh 1 {ci)"'dxh s ( c s 



■ exp 



-(X,(D + Q H )X) C 
la 



A 

2tt_ 
. (6.10) 
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Repeating the previous manipulations we then arrive at the explicit expres- 
sion 



<ii{n 1 )..4{n r )it hl {c 1 )...it hs {c s ) >= y \ [dA] exp[— 5m — Trln(C -1 (0)(f'(-4))] 



J[dipdip] n exp[-S N ]i/j(n 1 )...'il;(n r ) J 

d 

exp[— S x ] A (a2a i p 2 ) - ' 



dx 



d6 



d 



dXhAci)'" dxh s ( c s) 



2tt 
(A-X) . 



(6.11) 



7 Quark confinement 



Let us consider the quark propagator 

< X(x)T(x,y)X(y) >=\\ [dU] J [dA] exp[-S YM - S M ] 
J [dXdX] X(x)T(x,y)X(y) exp[-S N - S c - S Q ], 



(7.1) 



where T(x, y) is an arbitrary string of gluons making the above integral gauge 
invariant. We can repeat the previous manipulations until we arrive at 



< 



X(x)r(x,y)X(y) >=\\ [dA] ex V [-S M - Tr ln(C'- 1 (0)C(A))] J [dXdX] r 



expf-^A? - S Qn ] J 
J [dXdX] c X(x)T(x, y)X{y) exp 



dx 



dQ 



exp[-S x ] / [dU] exp[-S Y M] A 
2ttJ c J 

1 _ 

-(A, (D + Qh)X) c — SQ nc — Sq c 

Gj 



(7.2) 



In the above formulation the quarks appear in interaction with the auxiliary 
fields and the nucleon fields, and have an effective mass equal to k n (p 2 <f) + 
\/VLm) ~ \fVLpk n / y / a2a t . This mass seems to diverge in the continuum limit, 
but we will see that the constant k n must be tuned with the lattice spacing in 
a way which will not be determined in the present paper. What is relevant, 
however, as far as the perturbative expansion is concerned, is the dependence 
on Q. Since the quark effective mass in the chirally broken vacuum grows 
like x/T), the saddle point expansion results to be a strong coupling expansion 
for the quarks, so that if the sites x, y are n lattice spacings apart from one 
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another, the first nonvanishing contribution to the quark propagator occurs 
at an order of the expansion not smaller than n. As a consequence the quarks 
are never produced to any finite order, a fact first observed in [7|] . In this sense 
we can say that they are confined. Since the quark effective mass to leading 
order in Q does not depend on the breaking parameter, this confinement is a 
genuine consequence of the spontaneous chiral symmetry breaking, in whose 
absence there would be no saddle point expansion. 

In conclusion in the present approach the quarks do not propagate at 
the perturbative level, so that they might become alive only due to nonper- 
turbative effects. The situation is reversed wr to the standard perturbation 
theory where the quarks appear as physical particles whose confinement is 
intrinsically nonperturbative. 

Since the quarks have no poles whatsoever, it seems that we do not have 
to worry about the spurious ones. If this turns out to be confirmed by the 
study of the anomaly in the present approach, we can altogether forget the 
Wilson term for the quarks. In any case we can safely assume the Wilson 
parameter r g of order 1 / and neglect the quark Wilson term to leading 
order in our expansion. 



8 The nucleon-pion interaction 

As an application we investigate the pion-nucleon interaction to lowest order. 
It occurs in Xi j2 ,o , the lowest order term of Zi j2 . Since the integrand will 
result to be independent on the gluon field, the corresponding integral will 
factor out and we will therefore omit it from the start 



Zi,2,o = o(*SQn) 2 A / [d\d\} c SQ nc S Qcn exp 



-(A, DX) ( 



-a 3 a t (SQ T 



c,e,f] 



{c-e et )X hl (c + e, 



(8.1) 



We look for the terms in Ii,2,o proportional to ip, ip, namely containing three 
A at one and the same site and also three A at one and the same, possibly 
different, site. They are obtained by retaining only the contributions with 
e = — rj, and only the mass term from Sq u 

3-1,2,0 — 7" "/ "'(.) / , \ Vr + , LJ <■ Vr. 



-^^WqY. (Qc+e et ,cD c X Q C 



hjhi 
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(%x A/i 2 \h 3 A/n Afc 2 \ h3 ) c+eet 



h\h\ 



hh lh2 h 3 fi h^hMc + e et ) V/(c + e et ). (8.2) 

We need the transposed of the product QD~ X Q. Assuming r q of the order 
1 / y/Ti we can neglect the Wilson term of the quarks so that 

T 11 

{Qc+eucD^Qcc+e) = jj-^CT 2lt {v t ) + (c)D c 

lt v t (c)C- X T 2 <g> I c , (8.3) 
where I c is the unit in color space. Now by using the identities 



h h! h 2 h 3 flhih 2 h z (vf ) r i {V+ ) T2 (V t + ) T3 V tTl V tT2 V tT:j 



where 
we get 



flh^hsfihihzhsiV^hVtr, no summation (8.4) 
V tT (c) = exp(ie T a t A t (c)), (8.5) 



3 - 1 

2i >2>0 = -a^mJS^+^^cTrfC- 1 }^ 



^Jhifehs ^hift2hs [V'(c + e rf )V + (c)]_[l/(c)^(c + e et )] / . (8.6) 

The sum over colors has been performed. Using the explicit expressions of 
the functions h 

3 1 1 



8 24 2 -'"-V"" ^£> CJ D+ 

^(c + e et )U + (c) [3L> lc + i£> 2c ] V(c)^(c + e et ), (8.7) 

where we have decomposed D according to 

D = D 1 +iD 2 . (8.8) 
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Retaining the leading terms in in the matrix D and summing over e we 
get the final result 



Zi,2,o = ig N ~ w a 3 2a t J2^( c + e t)V + (c)>y 5 T ■ x(c)V(c)^(c + e t ) 

c 

+a 3 2a t 5m N ifj(c + e t )ip{c + e t ), 



(8.9) 



where 

3 1 2a t a 2 m 2 Q 
is the nucleon-pion coupling constant and 



5m 



N 



9 1 2a t arriQ 



~ m Q 



(8.10) 



(8.11) 



8 24 2 V a y/Upk N K 

is a renormalization of the nucleon mass. In the continuum limit the shift in 
the arguments of the fields can be ignored, and the link variables V disappear. 

The interpretation of Qn-tt as the pion-nucleon coupling constant follows 
by the evaluation of the 3-point correlation function 

< ^(mMn 2 K(c) >=\\ [dA] exp[-S M - Tr ln(C'- 1 (0)C(A))] 



exp[— S x ] A (a2a t p 



exp[-S , A r]^(ni)^(n 2 ) 
d 



2\-l 



dXh{c) 



2tt 

(A- 1 Xi,2,o) 



de 



(8.12) 



Since the term with the derivative of 2i,2,o does not contribute in the contin- 
uum limit, we have 

< ?(n!Mn 2 H(c) >=\\ [dA] exp[-S M - Tr ln(C- 1 (0)C(A))] J [##] n 
i){ni)i){n 2 ) exp[-SW] J 



Xh(c)T 1>2 , exp[-S x ]. 



(8.13) 



9 Conclusions 

We have derived a perturbative expansion in QCD in terms of the quark 
composites with the quantum numbers of the nucleons and the auxiliary 
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fields which replace the chiral composites. The coefficients appearing in the 
perturbative series are defined in terms of integrals over the quark and gauge 
fields. 

The expansion parameters are the inverse of Q, the number of quark 
components, the inverse of the constants k N ,k n , entering the definition of 
the composites and the quark mass tuq. The nature of the expansion is 
however at present not fully understood: k N is indeed accompanied by large 
numerical factors related to the number of quark components, but we do not 
know how the powers of these factors will be related to the powers of A; at in 
higher order terms. We do not see at present any other way of clarifying this 
issue than by studying higher order terms. 

To leading order the nucleon-pion interaction does not involve the gluon 
field, so that it is not affected by the gluon dynamics, but it is solely deter- 
mined by the hadronic quantum numbers. We regard this as the result of the 
choice of the "right" variables, aimed at suppressing the degrees of freedom 
irrelevant in the regime of interest. 

One of the problems where the present formulation of QCD seems par- 
ticularly promising is the phase transition at high barion density ||. The 
reason is that in the usual way of introducing the chemical potential this 
is coupled to the quarks to allow an analytical integration to be performed 
on the latters. But the price to be paid is an interplay between gluon field 
and chemical potential which would not be there if we coupled the chemical 
potential to the nucleons, what we can easily do in the present formalism. 
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